NP563-82-O0dCR 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


CONTRACT  REPORT 


AN  ANALYSIS  OF  THE  WYNGAARD-LEMONE  MODEL  OF 
REFRACTIVE  INDEX  AND  MICRQMETEOROLOGICAL  STRUCTURE 
FUNCTIONS  AT  THE  TOP  OF  A  TURBULENT  MIXED  LAYER 

C.  W.  Fairall 

BDM  Corporation 

P.O.  3ox  2019 

Monterey,   CA     93940 

October   1982 
Report  for  the  period   15  April  -  30   Sept   1982 


Approved   for  public   release;    distribution  unlimited 


Prepared   for 


Naval  Environmental  Prediction  Research  Facility 
Monterey,   California     93940 


PedDocs 
D   208.14/2 

NPS-63-82-006CR 


f       0«-  -8a-  oobc&. 


IM.V-  &»-  <** 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California  93940 


D.  A.  Schrady 
Rear  Admiral  J.  J.  Ekelund  Provost 

Superintendent 

Tae  work  reported  herein  ~2~»*~J*£i  SS^^^t 
Prediction  Research  Facility  by  BDM  &^^  *^L^*^  "Optical  tur- 
Number  N00014-a2-C-0251.  The  work  was  part  of  a  P£*^«££ onmental  Pre- 

s:  s^rr^r^w5^^^*^^.  *.  l.  «* 

Reproduction  of  all  or  part  of  this  report  is  authorized. 
This  report  was  prepared  by: 


Unclassified 


KHS-SWCS  VO  A3 

m  ioohos  3ivnavaoisod  ivavn 
Afcvysn  xonx  A3iana 


SCCUmTY   CLASSIFICATION   Of   THIS  PACE  'W>«n  Oa(«  £nf*r«d) 


REPORT  DOCUMENTATION  PAGE 

I.     REPORT   NUMBER 

NPS63-32-0Q6  CR 


— 


READ  INSTRUCTIONS 
BEFORE  COMPLETWr,  KORM 


2.  GOVT  ACCESSION  HO.     1.     RECIPIENT'S  CATALOG  NUMBER 


5.    TYPE  OP  REPORT  ft  PERIOD  COVEREO 

Contractor  Report 

15  Apr  82  -  31 Se?  1982 


TITLE  (■•«<*  5uo.ii/.) 

An  Analysis  of  the  Wyngaard-LeMone  Model  of 
Refractive  Index  and  Micrometeorological  Structur 
Functions  at  the  Top  of  a  Turbulent  Mixed  Layer 


*.    PERFORMING  ORG.   REPORT  NUMBER 


7.    AUTHORS  •> 

C.    W.    Fairall 


»      CONTRACT  OR  GRANT  NUMBER!**; 

JT00014-82-C-O251 


>.     PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 


BUM  Corporation 
P.O.  Box  2019 
Monterey,  CA  93  940 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

Work  Order   1202 


11.     CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval  Environmental  Prediction  Research  Facility 
Monterey,  California  93-940. 


12.     REPORT  DATE 

October   198: 


13.     NUMBER  OF  PAGES 

101 


14.    MONITORING  AGENCY  NAME  I   AOORESS</f  OllUtmnt  ttam  Contwoillnm  Otllem) 


IS.     SECURITY  CLASS,  (at  rma  rmpori) 

Unclassified 


ISA     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


14-     DISTRIBUTION  STATEMENT  (at  tttim  Assort) 

Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (of  tft.  imtrKi  mntoroo)  In  Block  20,  II  <Httormn  /ram  Rmport) 


I*.     SUPPLEMENTARY  NOTES 

Principal  Investigator,  K.  L.  Davidson,  Professor  of  Meteorology, 
Naval  Postgraduate  School,  Monterey,  CA   93940 


19.     KEY  WOROS    Cormnua  an  rartmrn  «<a»  .•/  noemmmmrr  *id  lamntltr  *T  3/ocjr  manOmr) 


20.     ABSTRACT  (Canitnum  an  rwrmtmm  too  It  imchmtt  and  lamnttty  ay  oioclt  .iuoo«r) 

"  The  Wyngaard  and  leMone  C198Q)  model  of  interracial  turbulence  structure 
functions  (temperature,  C<j  ,  and  water  vapor,  Cq  }  in  an  entraining  mixed 
laver  is  analyzed.   The  model  indicates  that  in  the  interracial  region  Z    as  ZJ ) 
C_  is  proportional  to  (£X\     Z.    -'  •>   9  4/£G  where  X  =  T  or  0,  IX  is  the 


is  the  neight  of  zhe    interface,  and  9_rj(c 


jump  in  X  across  the  interface,  Z 

the  convective  mixed  layer  scaling  parameter  for  temperature.   Although 

based  on  a  number  of  assumptions  (referred  to  as  the  "quasi-steady" 


00 


FORM 

JAN  73 


1473 


EDITION  OF    t  NOV  6S  IS  OBSOLETE 
S/N    U02-014-4601 


Unclassified 


SECURITY  CLASSIFICATION  OF   THIS  PAGE  (Whmn   Dm,,  inf 


Unclassified 


-LLUWITY  CLASSIFICATION  OF  THIS  PAGEfWh**  Dmta  Enttrmd) 


approximation),  the.  model  is  found  to  have  more  general  application.   A 
theoretical  analysis  indicated  that  the  model  might  not  apply  where  A0 
is  large  (on  the  order  of  10  Kl ,  particularly  for  C  z.     A  comparison 
against  23  aircraft  profile  measurements  revealed   that  the  model  agreed 
within  a  factor  of  three. 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  f»AGE(TW««n  Oatm  ZnfMl 


ABSTRACT 

The  Wyngaard  and  LeMone  (1980)  model  of  inter facial 
turbulence  structure  functions  (temperature,  Gp2#  and  water 
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model  indicates  that  in  the  inter  facial  region  (Z  ~  Zj_)  Cx2 
is  proportional  to  (AX)2  Zj_"2/3  9V*/A9V  where 

X  =  T  or  Q,    AX  is  the  jump  in  X  across  the  interface,  Zj_  is  the 
height  of  the  interface,  and  9V*  is  the  convective  mixed-layer 
scaling  parameter  for  temperature.   Although  based  on  a  number  of 
assumptions  (referred  to  as  the  "quasi-steady"  approximation),  the 
model  is  found  to  have  more  general  application.   A  theoretical 
analysis  indicated  that  the  model  might  not  apply  where  A0V  is 
large  (on  the  order  of  10  K) ,  particularly  for  CT2 ,   ^ 
comparison  against  23  aircraft  profile  measurements  revealed  that 
the  model  aareed  within  a  factor  of  three. 
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I   INTRODUCTION 

This  report,  is  a  theoretical  and  experimental  analysis  of  a 
model  (Wyngaard  and  LeMone,  1980)  used  to  calculate  the  refractive 
index  structure  function  parameter,  Cn2 ,  at.  the  inter  facial 
region  at  the  top  of  an  entraining,  turbulent  mixed  layer. 
Cn2  is  related  to  the  micrometeoro logical  structure  functions 
for  temperature,  C-j-2  #  humidity,  Cq2,  and  T-Q  covariance, 
C-J.Q.   The  mixed  layer  inter  facial  region  is  important  for  EM 
propagation  because  Cn2  is  greatly  increased  by  large  T  and  Q 
fluctuations  due  to  the  entramment  of  warm,  dry  air  from  the 
nonturbulent  atmosphere  above  the  mixed  layer. 

Assuming  that,  the  rate  of  entrainment  is  in  equilibrium  with 
the  free  tropospheric  virtual  potential  temperature  (buoyancy) 
lapse  rate,  the  model  indicates  that  C^   is  proportional  to 
(AX)2  9V*  Zi-2/3/£0v  where  X  is  T  or  Q,  AX  the  jump 
at  the  interface,  Zj_  the  height  of  the  ooundary  Layer  and 
9V*  the  convective  temperature  mixed  layer  scaling  parameter. 
The  theoretical  basis  of  this  model  is  examined  and  four  data  sets 
from  the  MPS  aircraft,  measurements  program  are  used  to  t.est.  the 
mode  1. 


II   THEORY 

A.   Background 

The  structure  function  parameters  for  temperature,  C«p2# 
and  specific  humidity,  cq2#  are  ro  be  evaluated  in  the 
inversion  region  by  averaging  between  heights  Z   =  hQ  ana 
Z  =  h2  (see  Fig.  1).   The  complete  theory  was  developed  by 
Wyngaard  and  LeHone  (1980),  hereafter  referred  to  as  WL,  so  only  a 
summary  of  the  derivation  will  be  presented  in  this  report.   In  a 
few  instances  WL's  work  will  be  expanded  to  make  certain 
assumptions  and  manipulations  more  explicit. 

The  height  h0  is  defined  as  the  top  of  the  mixed  layer 
where  wdv  =  0.   At  1*2  both  fluxes  and  flux  divergences  are 
equal  to  zero.   The  average  structure  functions  are 

<CT2>  =  Ah-1   /   CT2  dZ  (la) 


''O 


<CQ2>  =  Ah-1  J 


Cq2  dZ  (lb) 


where  Ah  =  Y\2    -   hQ  and  the  0,  2  on  the  integral  denotes  h0,  h2  • 

The  average  structure  functions  are  related  to  their 
respective  dissipation  rates  by  r.he  Corrsin  equation 

<CT->  =  1.6<£>-l/3<Xg>  (2a) 

<Cq2>    =  1.6<s>"1/3<Xq>  (2b) 


where  £  is  the  rate  of  dissipation  of  turbulent  kinetic  energy,  X 
and  Xq  are  the  scalar  dissipation  rates  (the  factor  1.6  implies 
X   is  the  rate  of  dissipation  of  temperature  variance  3^). 

B.   Evaluation  of  <X> 

For  the  moment,  the  development  will  be  confined  to  the 
specific  humidity  (Q) .   The  dissipation  rata  is  calculated  from 
the  scalar  variance  budget  equation  (Q  denotes  mean  while  q 
denotes  fluctuating  specific  humidity;  later  in  the  paper  q  will 
denote  mixing  ratio,  Q/p)  . 

dv/dt  +  Wdv/dZ  +  d(wq)/dZ  +  2  owq  d(Q/p)/dZ  =  -Xq         (3) 

where  v  =  q2,  w  is  the  mean  vertical  velocity  (subsidence)  and  p 
is  the  density  of  air.   Integrating  this  equation  from  hQ  to 
h2»  as  in  Eq .  1,  yields 

<Xq>  =  -<DQ>  -  <TQ>  -  <PQ>  (4) 

where  D  is  the  first  two  terms  in  Eq.  3,  T  (transport)  the  third 
and  ?  (gradient  production)  the  fourth.  Assuming  "quasi-steady" 
conditions,  >7L  show  that  <Dn>  and  <T,-i>  are  negligible  compared 


-         w 


r.o  <?->>:    therefore 


<  Xq  >   =   -  <  ?  q  >  (  5  ) 


At  this  point  the  generalized  inversion  structure  model 
(Deardorff,  1979)  is  introduced 

Q  =  Q0  +  Aq  f(z);  h0<Z<h2         (6a) 

dQ/dZ  =  AQ  df/dZ  (6b) 

where  f(Z)  describes  the  shape  of  the  Q  profile  in  the  inversion 
region  (assumed  to  be  the  same  for  G  and  T)  with  f(hQ)  =  0  and 
f(h2)  =  I,    Q0  is  the  mixed  Layer  value  and  Aq  the  jump  in  Q 
across  the  inversion.   Substituting  Eq.  6b  into  Eq.  5  and 
integrating  by  parts  one  obtains 

-<PQ>Ah  =  2  AQ   I   d(wq)/dZ  fdZ  <7) 

-'O 

The  mean  Q  continuity  equation 

-d(wq)/dZ  =  dQ/dt.  +  W  dQ/dZ  (8) 

is  used  in  2a .  7  to  obtain 


r2  rl 

-<?.->  Ah   =    -2_C     I     dQ/dt    fdZ    -    2  AC  'v    iQ/dZ    f   dZ  (9) 

Jo  Jo 


Th< 


dQ/dt   =    iQn/dt    +    f    dAQ/it  (10 


and  Eq.  6b  can  be  substituted  into  Eq.  la.   First  the 

"quasi-steady"  assumption  is  invoked,  setting  the  following 

conditions 

dAQ/dt   =  0  (11a) 

dA9v/dt  =  0  (lib) 

dhQ/dt   =  0  (He) 

dAh/dt   =  0  (lid) 

However,  since 

dhQ/dt  =  WQ  +  WeQ  (12) 

then  Eq .  lie  implies  Weo  =  -WQ.   Assuming  constant 
divergence 

W  =  W0  Z/hQ  d3a) 

dW/dZ  =  w0/h0  (13b) 

w2  -  (l  +  a)  wQ  (13c) 

where  a  =  Ah/hQ  is  the  normalized  thickness  of  the  interracial 
region.   Employing  these  relations  in  Eq.  9  and  doing  the  second 
integral  by  parr.s  gives 

-<?~>Ah  =  -2  AG  hYQ  iG0/dt  +  (  AQ  )  2W9Q  (I  tci-oe^)  (14) 


where  the  interracial  functions  Yq  and  3,^  are 


Yq  =  Ah-1  Jo  f  dZ  (15a; 

-4 


o 

•2 
/o 


Zo  =  ^h"1  L  fz   dZ  (15b; 


The  time  derivative  term  in  Eq.  14  is  eliminated  by- 
integrating  the  conservation  equation  (Eq.  10)  from  hQ  to  V12 


(16 


AhdQQ/dt  -  Weo  AQ(l  +  a  -  ayQ )  =  wq0 
which  is  substituted  into  Eq.  14  to  obtain  (WL  Eq.  42) 
-<PQ>Ah  =  -2AQYQ  wo~  +  ( AQ)2Weo[-2YQ(l+a-aYQ)+(l+a-azQ)]      <17) 
Later  in  their  paper,  WL  use  the  equation 


-AC  Weo  (1  +  a  -  aYQ)  =  wq0 


(13) 


which,  in  view  of  Eq.  14,  obviously  implies  dQ0/dt  =  0.   Since 

WL  have  already  required  that  dAQ/dt  =  0,  this  solution  appears  to 

be  quite  restrictive.   If  Eq.  18  is  used  in  Eq.  17  then 

-<?Q>Ah  =  (AC)2  Weo  (I  +  o  -  2Zq)  <19> 

Despite  the  simplicity  of  Eq .  19,  WL  prefer  to  keep  the  vq0 
term  separate  in  r.heir  development.   The  primary  reason  for  this 
is  to  simplify  the  analagcus  development  for  9V  since 
wyvo  =  0.   Therefore,  WL  now  employ  r.'ne  "quasi-steady" 
entrainment  formula 


W   =  0.8  W*  S-l/(l+a)  (20) 

where 

S  =  gr92  h2/(w*2T)  (21) 

with  Tq2  =  d9v/dZ  at  Z  =  h2  and  W*  is  the  convective 
scaling  velocity  (Zj_  =  ho/0.8) 


W*3  =  g  w6vs   Zi/T 


(22) 


and  "s"  denotes  the  surface  value. 

Rather  than  make  an  explicit  substitution  for  Weo  at  this 
point,  one  could  keep  Weo  as  a  variable,  giving 

<XQ>  =  -2AQ  YQ  wo~/(ah0)  +  (AQ)2  Weo (1 +a)SFQ/h0  (23) 

where 

FQ  =  C-2YQ  (1+a-aYg)  +  ( l+a-aZQ)  ]/  (  a  ( l+a)S  )  (24) 


Using  the  WL  solutions  obtained  for  Eq.  24  (and  Yq  =  1/2) 

FQ    =     (6R)"1 


(25) 


26) 


then 

<Xq>    =    -AQ   wq*c    YQ/(ah0)    -    (AC) 2Weo(l+a)S/ (6Rhc) 

Where  R  -   g    A0V  hQ/(W*2T)  (27) 

The  final  result  is  obtained  by  substituting  for  vq"c  in  Sq. 
26.   First,  the  9V  continuity  equations  ar  "n0  and  ar.  h2  are 


combined  with    the  h0   to   h2    integral    form   similar   to    Eq.    16   to 
oroduce   the   relation 


d    A<Vdt   +   '*eo    Ah(l+a)r92    -  Weo    AGv(l+a-aYQ)    =  w6vo  (28) 

Since   WL   assume   dA0v/dt   =   w6vo    =    0, 

A0v(l+a-aYQ)    =   "92      h(l+a)  (29) 

Using   Eq.    29    and   Eq .    18    in    Eq.    26,    one   obtains 

<XQ>    =    (Aq)2   weo(l+a)s(l+6-1)/(h0R)  (30) 

Mote  that  the  first  term  in  Eq.  26  (which  was  proportional  to 
wqQ )  is  six  times  as  large  as  the  second  term  (WL  obtain  15/2 
for  this  ratio  because  they  use  two  separate  formulae  for  Weo 
which  differ  by  a  factor  of  4/5,  i.e.  6*5/4  =  15/2). 

The  development  for  temperature  is  parallel  until  the  9V 
aquation  analagous  to  Eq.  26  is  reached.  Since  w6vo  =  0,  the 
final  result  is 

<X9>  =  (A3v)2(i+a)weo  S/(6Rh0)  (31a> 

<Xq>  =  7(AQ)2(i+a)wso  S/(oRhc)  <31b> 

C.   Structure  Functions 

The  final  step  in  r.his  process  is  ro  specify  that  <  s>  is  one 


half  the  value  typically  found  at.  Zj_  under  convective  conditions 
<£>l/3  =  (0.2)V3  w*  Zt-l/3  (32) 

Assuming  the  "quasi-steady"  entrainment  rate,  the  structure 
functions  become 

<C02>  =  3.9(AQ)2  8v#/(Zi2/3^0v)  (33a) 

<CTv2>  =  O.5A0V  9V*/Zi2/3  (33b) 

where  8V*  =  w6~vs/W*.   The  virtual  temperature  structure 
function  is  related  to  the  temperature  structure  function, 
CT2,  by 

<CT2>  =  2  Ti    <CT2>/A0v  (34) 

where  T^  is  the  function  given  by  WL. 

One  point  worth  more  discussion  is  the  approximation 

Fq  =  (6R)"1  and  the  final  forms  of  Eq.  31.  Suppose  the 

results  of  Sq.  19  were  used  and  a  different  function  defined 

<Xq>  =  (AQ)2WQO(l+a)  5  GQ/hQ  (35a) 

<X9>  =  U0v)2Weo(l+a)  3  FQ/h0  (35b) 

where  Fo  remains  as  oer  Sq.  19  but. 


GQ  =  (1-kx-  azQ)/(a  (l+a  )s)  (36) 

Using  Eq.  29  one  can  show 

Fq  =  GG  "  2yq/r  (37) 

Following  the  calculations  WL  have  in  their  Appendix  A,  a , 
Fq  and  Gq  are  unique  function  of  R/'S  (providing  the  mixed 
layer  gradient  is  zero) .   Gq  is  considerably  less  variable  than 
Fq.   The  following  formula  are  reasonable  approximations  for 
0.1<R/S<10 

a  =  0.96  R/3  -  0.11  (R/S)1*5  (38a) 


Gr  =  (1  +  0.064  '/R/SjR"1 


(38b; 


FQ  «  (1  +  0.28  n/r/S)R-1/6  (38c) 

These  formulae  lead  zo   slight  modifications  to  the  structure 
function  equations 

<CQ.2>  =  3.3  (AQ)2ev*/(Zi2/3A0v)  (39a) 

<CT.2>  =  Q.57  (A9V)9V*DT/Zi2/3  (39b) 

where  Om  =  1  +0.22  yR/S. 


The  equations  for  Cq2  anc}  q^2   can  <oe  written  in 
various  general  forms  (now  dropping  the  bracket  notation) 

CX2/((AX)2DXEX)  =  1.14  9v*/(AGvZi2/3)  (40) 

or,  without  substituting  explicitly  for  Weo  and  e 


->      o        0.53(l+a)  rw  , ,, , 

CX2/((AX)2DXEx)  =  A(j   <£>i/3?  (41) 


where  Dq  =1,  Eq  =  3,  and 

ET  =  Ti/A0v  (42) 

D.   Discussion 

It  is  of  interest  to  ponder  the  significance  of  the  various 
"quasi-steady"  assumptions  (Eq.  11,  12,  13).   Suppose  we  exhume 
the  original  conservation  equation  integrals  from  Deardorff's  (1979) 
paper  (his  Eq .  13  and  21).   Assuming  only  horizontal  homogeneity 
and  constant  divergence,  the  general  aquations  become 

Ah  3QQ/dt  +  AhZQ  iAQ/dt  -  1QL  (  I  -  YQ)  Wa2  +  -Qweo^  =  ~-c  <43a) 

Ah  ievo/it  -  AhYg  dA0v/dt  -  A0V[(1  -  YQ)  We2  +  *q  ;^o^  =  °    <43b> 
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Thus  Eq  16  and  the  9V  analogue  can  be  reproduced  by  requiring 

We2  =  (1  +  a)  Weo  (44a) 

dAQ/dt.  =  dA0v/dt  =  0  (44b) 

It  is  not  necessary  to  require  We  =  -W,  dhQ/dt  =  dh2/dt  =dAh/dt  =  0 
This  explains  why  WL  found  excellent  agreement  with  Aschurch  data 
where  W  =  0  and  Weo  ~    10  cm/s. 

Similarly,  the  general  forms  for  the  dissipation  integrals 
are 

-<PQ>Ah  =  -2 Ah  AQYQdQ0/dt  -  2AhZQAQdAQ/dt 

+(AQ)2[(i  -  Zq)  We2  +  ZqWso]  (45a) 

-<?9>Ah  =  2AhA0v  (Yq2  -  Zq)  dAGv  /dt 

+  (A0V)2[-2YQL(1  -  YQ)  We2  +  YQ  Weo]  (45b) 

+  (1  -  Zq)  We2  +  Zq  Weo] 

which  reduce  to  Eq  14  if  rhe  conditions  of  Sq  44  are  met.. 

Since  entrainment.  and  surface  flux  tend  ro  counteract  each 
ether  in  the  Q  case,  it  seems  quite  reasonaoie  to  assume  -hat  the 
iAQ/dt  and  dQ0/dt  -arms  are  negligible  in  Eq  45a 

-<?Q>Ah  =  (1Q)2  [(I  -  ZQ)  We2  +  Zq  WeQ]  (46) 

Instead  of  making  the  assumption  Sq  44a,  suppose  ve  simply  assume 
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we2  =  weo  =  a'^vs/A0v 


(47) 


which  is  the  standard  cloud-free  result  from  Lilly  (1968)  where 
typically  a  =0.2.   Then  one  can  easily  show  that 


<CQ2>  =  3.3  (AQ)2  9v*/(Zi2/3A0v) 


(48) 


which  is   identical  to  the  WL  result  as  expressed  in  Eq.  39ai   In 

other  words,  the  combination  of  "quasi-steady"  assumptions 

We2  -  (1  +  a )Weo  and  We2  =  w6vs/( r9hQ)  are 

equivalent  to  the  assumptions  of  Eq.  47  even  though  they  may  imply 

vastly  different  entrainment  rates . 

If  one  uses  the  assumptions  of  Eq.  47  and  parallels  the  WL 
development,  then  the  equivalent  to  Eq.  18  is 

wqQ  =  -  AQ  Weo  (49) 

and  the  equilibrium  condition  form  the  3V  equation  is 

«6vo  =  weo  (ahor9  ~  -°v)  (50) 

which,  assuming  w8 vo  =  0,  gives 

a  =  R/S  (3D 

The  results  for  9V  are  also  interesting  cecause  it  is  not 
clear  that  the  dAGv/dt  term  should  "ce  negligible  compared  to  the 

1  "3 


other    terms    in   Sq   46   b.       Suppose   we    Let 

-<Pg>  Ah    =   A   +    3  (52) 

Then  the  dA0v/dt  term  is  small  if  A/B  is  small  (returning  to  the 
"quasi -steady"  format) 


hQ(YQ2  -  ZQ)  dA0v/dt  6(R/S)  (53) 

A./  B  = 


A0V  ( l  +  a)  w 


eo 


Since  Yq2  -  Zq  z    -0.1,  we  can  write 


-0.6  hn(R/S)   dA0v/dt 

A/B  =  2 Y. (54) 

(1  +  a)  weo  A0V 


The  magnitude  of  A/3  can  be  examined  by  using  the  general 
relation 

dA0v/dt  =  -d0o/dt  +  r9We2  (55) 

and  writing  a  simple  antrammenr.  formula  (e.g.  "quasi-steady"} 


*e2  =  '^vs-  (-  rqh0) 
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'he  integral  of  the  conservation  aquation  from  Z  -  0  to  Z  =  h0  gives 


h0  d0vo/dt  =  W9VS  +  weoAGv  (57) 

therefore 

dtQv   /dt  =  -  WeoASv/hQ  (58) 

using  Eq.  54  we  find 

0.6 

A/B  = R/S  (59) 

(1+a) 

A  good  example  is  the  Aschurch  data  quoted  by  WL  where  Eq.  57 
was  shown  to  be  applicable.   Since  R/S  =0.3  for  that  data,  A/3  = 
0.15  and  dA0v/dt  is  negligible. 

Certainly  the  conditions  set  by  WL  are  consistent  with 
neglecting  dA0v/dt.   It  is  not  clear  how  to  ident.ify  conditions 
where  this  assumption  is  invalid.   Eq .  54  cannot  provide  much 
guidance  because  it  is  based  on  solutions  to  Eq.  28  with 
dA0v/dt  =0.   It  is  interesting  that  in  the  conditions  where  the 
WL  equations  for  "quasi-steady"   entrainment  are  expected  to 
breakdown  (A0V  large,  R/S  >1)  then  r.he  Lilly  type  relations  give 
the  same  resulr.s  for  Cq2  .   if  the  dA0v/dt  terra  becomes 
important,  then  one  anticipar.es  the  WL  formulation  will 
underesr.imare  Gm^. 


Ill   ATMOSPHERIC  DATA 

A.   Measurement  Techniques 

The  measurements  were  made  using  a  single  engine  BelLanca 
Viking  aircraft  operated  by  Airborne  Research  Associates  of 
Weston,  MA.   The  instrumentation  and  data  processing  have  been 
previously  described  in  detail  (Fairall  et.  al.,  1980;  Schacher 
et.  al.,  1980)  so  only  a  brief  summary  is  given  here. 

i)     Mean  temperature,  T:   platinum  resistance  sensor  with 

standard  aircraft  mount. 

ii)    Mean  humidity,  Q:   cooled  mirror  dew  cell. 

iii)   Mean  windspeed,  U:   estimated  at  the  surface  from  the 

sea  state  and  DMV  navigational  aid.   The  present  LORAN  system 

was  not  available. 

iv)    Sea  surface  temperature,  Ts:   Barnes  PRT-5  IR 

radiometer. 

v)     C<t>2.   microthermal  sensors  (4.5  um  dia.  tungsten) 

in  the  paired  configuration. 

vi)    Cq^ :   Lyman-alpha  fast  humidiometers  using  the 

mertial  subrange  filter  method.   Absolute  calibration  based 

on  comparison  with  a  microwave  refractometer . 

vi i)   £:   hot  wire  (4.5  um  dia.  tungsten)  constant 

temperature  anemometer.   The  inertial  subrange  filter  method 

was  used. 
3.   Surface  Fluxes  and  Turbulence  Scaling  Parameters 

Surface  fluxes  were  evaluated  from  aircraft  measurements 
using  two  methods:  a)  bulk  aerodynamic  and  b)  dissipation 
(inertial  subrange).   The  fluxes  are  defined  in  terms  of  rhe 
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following  scaling  parameters: 


momentum:  Puws  =  -ou*^  (60a] 


sensible  heat:      pCp  w9s  =  -pCp  u*T*  (60b) 

moisture:  O^s  =  -pu*q*  (60c) 

The  momentum  flux  is  also  referred  to  as  the  Reynolds  stress,  t. 
Note:  the  bulk  method  was  not  used  overland. 

1.   Bulk  aerodynamic  method. 

The  exact  details  were  described  in  a  recent  paper  (Davidson 
et  al,  1981).   Using  Eq.  4a  from  that  paper,  one  can  relate  the 
values  of  some  meteorological  variable  (temperature,  moisture  or 
wind  speed)  at  the  sea  surfce,  Xs,  and  at  some  height  Z  in  the 
surface  layer,  X2,  to  the  scaling  parameter,  X*: 

u*  =  uzk[2n  (Z/Z0)  -  ?u  (Z/L)]"1  (6Ia) 

T*  =  (Tz  -  Ts)  SkC  Zn  (Z/ZoT)  -  >?m  (Z/L)]"1  (61b) 

q*  =  (qz  -  qg)  ?k  [in  (Z/ZoT)  -  ^T  (Z/L}:*1  (61c) 

where  ZQ  and  ZQ^  are  roughness  Lengths,  L  is  the  Monin-Cbukhov 
length,  3  and  k  are  constants,  and  ?u  and  %   are  empirical 
functions . 


2.   Dissipation  method. 

The  dissipation  method  relies  on  semi-empiricai  relationships 

of  inertial  subrange  turbulence  to  surface-layer  scaling 

parameters  (Fairall  et  al.  ,  1980).  The  equations  are 

u*  =  [(  e  k  Z)/$  (Z/L)]1/^  (62a) 

T*  =  [Z2/3  CT2/f(z/L)]i/2  (62b) 

Q*  =  [Z2/3  Cq2/(A  f(Z/L)]V2  (62c) 

where  £  is  the  dissipation  rate,  $  and  f  are  empirical  functions, 
and  A.  is  a  constant.   3ecause  the  structure- function  parameters 
Crp2  anj  Cq2  are  related  to  the  square  of  the  scaling 
parameter,  a  sign  ambiguity  exists.   This  can  usually  be 
eliminated  by  examining  the  low-level  height  dependence  of  £, 
Cq2  and  Crr.2  because  the  functions  $  and  f  have 
characteristic  profiles  for  stable  and  unstable  conditions. 

3oth  methods  yield  accuracies  on  the  order  of  10%  for  u*, 
jr0.02°C  for  T*  and  +0.02  g/m3  for  Q*  (note:  q*  =  Q*/p )  . 
G.   Data  Sets 

The  data  given  in  this  report  were  cctamed  in  four  field 
programs: 

i)    Panama  City  (?C),  1973  (more  detail  avaiiaole  in 
Fairall,  1979)  over  the  Gulf  of  Mexico  in  Florida. 

ii)   White  Sands  (WS ) ,  1979.   Two  profiles  over  the  desert 
under  highly  ccnvective  daytime  conditions. 


IS 


iii)  MA  GAT  (MG),  1980  (more  det.ail  available  in  Fairall, 
1980)  in  the  Monterey  Bay  area. 

iv)   Bahamas  (BH),  1980.   A  series  of  profiles  taken  near 
Andros  Island  in  the  late  fall. 

The  complete  data  sets  were  examined  to  remove  profiles  that 
encountered  boundary-layer  clouds.   A  total  of  23  profiles  were 
selected.   Graphs  of  the  mean  and  turbulence  profiles  for  each 
case  are  given  in  Appendix  A.   A  summary  of  the  basic  scaling 
parameters  is  given  in  Table  1. 
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TABLE  1. 


Meteorological  data  and  surface  scaling  parameters 
from  the  cloud  free  NPS  data  sets. 
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Site 

Date 

Time 

u* 

T. 

<?* 

1, 

A0v 

AQ 

a 

rfl 

-i 

ms 

K 

aka"1 

km 

K 

am" 

Kkm'1 

[_          J_ 

1               1 

1 

PC 

11/26 

1252 

.40 

-.082 

-.16 

.85| 

1 

-6.5 

.4 

5.5 

1      2 

PC 

11/26 

1436 

.23 

-.095 

-.16 

.90| 

.5 

-2.3 

.1 

5.3 

1      3 

PC 

12/2 

1405 

.24 

-.  14 

-.18 

.231 

4 

-.5 

.7 

4.6 

1      4 

PC 

12/10 

1324 

.33 

-.35 

0 

•  9l| 

6 

-1 

.35 

10 

|      5 

PC 

12/10 

1410 

.32 

-.49 

-.49 

•75| 

.3 

-1 

.15 

11 

6 

PC 

12/10 

1523 

.34 

-.48 

-.48 

.85 

3 

-1.3 

.25 

11 

7 

PC 

12/10 

1637 

.34 

-.49 

-.50 

1 

.1 

3 

-3 

.1 

17.5 

1      8 

PC 

12/11 

1021 

.28 

-.44 

-.43 

.7 

3 

-1 

.5 

9.5 

9 

PC 

12/13 

1154 

.19 

0.21 

-.47 

.6 

1.5 

.2 

.35 

10 

10 

PC 

12/13 

1459 

.17 

-.20 

-.42 

.5 

.5 

-2 

.4 

11 

11 

ws 

10/17 

1330 

.47 

-.42 

0 

1 

•1   1 

1.5 

-2.5 

.1 

3.0 

12 

ws 

10/18 

1330 

.47 
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0 

1 

.9 

1.5 

-2.5 

.1 

3.3 

13 

MG 

4/30 

1610 

.23 
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-.11 

.  36| 

6.5 

-4.5 

.35 

9 

14 

MG 

5/4 
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.36| 

11 

-5.2 

.4 

10        1 
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MG 
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1201 
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-.075 
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.461 

9 
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.5 

15 

16 

MG 
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.54| 

9 

-5 

.2 

15 
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MG 
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.41 

-.04 
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.  23J 

7 

-2 

•5 

9 

13 

3H 
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•  5    1 

1 

-2.5 

•p 

5 
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3H 
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.33 
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.551 

0 

0 

? 

4.8  1 

!   20 

3K 
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.90) 

2.5 

-8.5 

1   c 

6.3 

i  21 

3K 
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1 

.5   ! 

3.5 

-9 

.  15 

5  .  5   1 

j  22 

3H 
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1347 

.20 

-.16 

-.26 

1 

.  5 

3.5 

-9 

.3 

5.5   | 

123 

3H 

12  <'15 

1637 

.14 

-.14 

-.25 

1 

.  1 

1 

-4.  5 

.4 

6.3   ! 

2C 


IV.   RESULTS 

A  summary  of  the  secondary  scaling  parameters  used  for  the 
NPS  data  set  is  given  in  Table  2.   Also  shown  in  Table  2  is  a 
comparison  of  the  measured  and  model  assumed  values  for  e  at  the 
inversion.   With  very  few  exceptions,  the  model  assumption  (Eq. 
32)  is  very  good.   The  entrainment  velocities  calculated  from  the 
"quasi-steady"  assumption  used  by  WL  (Eq.  20)  and  the  more 
conventional  parameterization  of  Lilly  (1968). 

WQO/W*  =0.2  ev*/A0v  (63) 


are  also  calculated. 

In  Table  3  are  the  measured  values  of  Cp2  ancj  Cq2  at 
the  inversion  plus  their  normalized  forms 


Ix  =  Zi2/3  Gk2/((Ax)2  dx  Px)  (64) 

taken    from   Eq.    40.       According    to   WL    (Eq.    26),    the    theoretical 
value    is 

Ic   =    1.14    9V*/A9V  (65) 

which  is  the  same  for  T  and  Q. 

A  direct  comparison  of  measured  and  calculated  values  of 
C.T.-2  and  Cq^  is  given  in  Fig.  2.   The  model  prauic-s  zhe 
measurements  within  a  factor  of  three.   The  uncertainty  is 
slightly  greaner  than  the  factor  of  two  suggested  by  WL  but 
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Table  2. 


Surface  scaling  (wevs  and  L ) ,  convective  scaling 
(W*,  9V*  and  £j_)  and  inversion  scaling  (R,  S 
and  Weo)  parameters.   Two  formulae  are  used  to 


estimate  W 
64. 


eo* 


"steady"  is  Eq  20  and  "Lilly"  is  Eq 
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Table  3. 

Measured  values  of  the  interfacial  structure  functions 
(CT2  and  Cq2)  and  their  resultant  values  for 

Ix  -  Z±2/2    CX2/((AX)2  DX  Fx) where  X^T  . 
or  U  Ttiese  are  compared  v/ith  theoretical  values,  Ic# 
using  the  "steady"  and  "Lilly"  entrainment  values. 
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includes  various  measurement  errors  and  uncertainties.   Mote  that, 
the  Cq2  data  has  a  greater  range  of  values  tnan  Op2.   This 
is  consistent  with  the  WL  model.   If  we  examine  the  function 


H  =  Zt2/3  CX2/(DX  9V*)  (66, 

then 

HT  =  FT  A0V  (67a) 

HQ  =  (AQ)2/A9V  (67b) 

A  graph  of  'drp  and  Hq  is  shown  in  Fig.  3  for  a  typical  range  of 
AG  v  anc*  AQ  from  the  MPS  data  set.  Note  that  Hrp  varies  roughly 
from  2  to  9  while  Hq  varies  from  4  to  72. 

The  entrainment  parameterization  was  tested  (Fig.  4)  by 
plotting  measured  values  of  Ix(Sq.  65)  against,  the  model  value 
( Eq .  66)  which  is  based  on  the  entrainment  formula  given  by  WL 
(Eq.  26).   This  plot  gives  a  much  higher  correlation  than  a 
similar  graph  (not  shown)  using  t.he  more  traditional  formula  due 
to  Lilly  (1968),  Eq.  62,  which  gives 

Ic«  (Lilly)  =  0.13  (l+a)r9  Zi   9V*/(A0V)2  (63) 

This  is  not  really  significant  because,  when  used  in  proper  com- 
bination with  Eq .  48,  the  Lilly  formulation  also  leads  to  Eq.  66. 

In  order  to  look  for  systematic  errors,  the  ratios  ( R^  and 
Rq )  of  measured  to  model  values  of  Z^~    and  Cq*   were 
calculated  and  plotted  against  AQv  (Fig.  5).   A  simple 
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log-average  yields  RT  =  1.15  and  Rq  =  1.3.   Figure  5a  weakly 
suggests  that,  the  model  underestimates  C^    (large  Rip )  when 
A0V  is  small  while  it  overestimates  when  A0V  is  large  (the 
Cq2  data  is  too  sparse  to  clear  up  this  question).   This  could 
be  due  to  an  error  in  the  estimation  of  A0V  and  AQ  (admittedly 
rather  subjective) .   An  examination  of  Fig.  3  suggests  that  a 
reasonable  adjustment  of  A0V  (several  tenths  K)  will  not  move 
the  data  points  substantially  closer  to  the  Rip  -  1  midline. 
Another  possibility  is  that  Eq.  20  tends  to  overestimate  WQO 
when  AGV  is  large  while  underestimating  for  small  A0V. 

Given  the  considerable  scatter  in  the  results,  the 
uncertainties  in  the  estimation  of  A0V  and  AQ  from  measured 
profiles  and  the  insensitivity  of  Cm2  to  A0V  and  AQ  it  is 
suggested  that  a  simplified  formula  for  Cm2  can  be  used  for 
application  to  radiosonde  quality  data.   If  one  assumes  (based  in 
Fig.  3)  that  EiT  -  5,  then 

CT2  ^  5  8V*  Z-l'2/3  (69) 

3ased  on  the  NPS  data  ser.  this  approximation  appears  to  be  at. 
leas-  as  accurate  as  the  more  complicar.ee  formula  (Fig.  6). 
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V   CONCLUSIONS 

The  Wyngaard-LeMone  inversion  Layer  scaling  has  been  examined 
theoretically  and  tested  against  a  data  set  obtained  by  NPS 
investigators  in  cooperation  with  Airborne  Research  Associates. 

The  theoretical  examination  indicated  the  following: 

i)    The  WL  theory  is  more  general  than  is  implied  by  the 
strict  assumptions  of  the  "quasi-steady"  theory. 

ii)   The  WL  development  can  oe  simplified  slightly,  leading 
to  modest  adjustments  of  the  normalization  constants. 

iii)  The  staady  state  assumption  that  dAQ/dt  is  negligible  is 
reasonable  under  most  conditions.   The  assumption  that  dAOv/dt 
is  negligible  may  not  be  justified  when  R/S  >  1. 

The  examination  of  the  atmospheric  data  indicated  the 
following: 

i)    The  assumption  that  z    at  the  inversion  is  proportional 
to  a  fixed  fraction  of  the  surface  buoyancy  flux  was  quite 
reasonable. 

ii)   The  WL  model  predicted  the  measured  value  of  CT2  and 
Cq2  to  within  a  factor  of  three. 

iii)  Some  evidence,  though  statistically  weak,  was  round  ^o 
suggest  the  model  overestimates  the  structure  functions  for  large 
A0V  (>  3K; )  while  it  underestimates  ior  snail  A0V  ( <  2K) .   Dn 
the  other  hand,  this  could  be  a  manifestation  of  the  Stein  affect. 
for  comparison  of  iar.a  sets  subject  ro  error  vhere  small  values 
are  usually  overestimated  and  large  quantities  are  usually 
underestimated . 


Based  on  these  results,  it  is  obvious  that  a  major  weakness 
of  the  model  is  its  reliance  on  an  entrainment  formulation  that  is 
too  restrictive.   The  two  extremes  of  the  buoyancy  jump  (A0V) 
may  involve  different,  entrainment  regimes  (e.g.  encroachment, 
convective  instability  or  the  Lilly  formulation).   It  would  also 
be  useful  to  include  the  effect  of  inversion  windshear  on  We  and 
on  the  structure  functions.   Another  area  of  investigation  might 
be  stable  surface  layers.   These  may  be  very  important  for  surface 
optical  propagation  because  C<p2  values  are  often  sizeade  and 
Z^  is  usually  small  (on  the  order  of  100m). 
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APPENDIX  A 

This  appendix  contains  graphs  of  mean  (6V,  q)  and 
turbulence  (C<p2/  Cq2,  e)    profiles  for  each  of  23  data 
sets.   The  site  designations  are  defined  in  Section  III-C.   The 
abstraction  of  this  data  to  obtain  the  relevant  parameters  (Tables 
1,  2,  3  in  the  main  text)  is  described  in  Sect.ion  III. 
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